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Abstract 

In a recent paper [1] we have proposed a new approach for extracting the wave function 
of the TT- meson (Ptt{x) and the masses and wave functions of its first resonances from 
the new QCD sum rules for non-diagonal correlators obtained in [2]. Here, we test our 
approach using an exactly solvable toy model as an illustrating example. We demonstrate 
the validity of the method and suggest a pure algebraic procedure for extracting the 
masses and wave functions relating to the case under investigation. We also explore the 
stability of the procedure under perturbations of the theoretical part of the sum rule. In 
application to the pion case, this results not only in the mass and wave function of the 
first resonance (vr'), but also in the estimation of 7r"-mass. 
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1 Introduction 



An important problem in the theory of strong interactions is to calculate, from the first prin- 
ciples of QCD, hadronic wave functions (WFs) (pn{x), ipn{xi,X2,X3), . . . [3] that accumulate 
necessary information about the non-perturbative long-distance dynamics of the theory. These 
phenomenological functions naturally appear as a result of applying "factorization theorems" to 
hard exclusive processes [3], [4], [5]. It seems now that only the QCD Sum Rule (SR) approach 
[6] and lattice calculations [7] can provide information about the form of hadronic WFs. 

The most popular set of hadronic wave functions, due to V. L. Chernyak and A. R. Zhitnits- 
ky (CZ) [8], was produced with the help of QCD SR for the first moments {C^)h {N = 2, 4) of 
mesons WFs. But the method they used has three main drawbacks: 

1. reducing the initially non-local objects to a few local ones: (g(O)g(O)), (G(O)G(O)), . . ., - 
leads to the fast growth with N of the power corrections in the operator product expansion 
(see discussion in [12]); 

2. the reconstruction of the whole function (pn{x) from just two non-trivial moments (C^)7r 
and is an unreliable procedure; 

3. the phenomenological model of the spectral density, based on local quark-hadron duality, 
is too crude for QCD SRs for the WF. 

Now it is known that the hadronic WFs are rather sensitive to the structure of the non- 
perturbative vacuum [9]. Therefore one should use a non-local condensate like (g(O)-E'(O, z)q{z)) 
and/or {G{0)E{0, z)G{z)) [10] which can reflect the complicated structure of the QCD vacuum. 
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The modified QCD SR with non-local condensates have been constructed in [9] and it has 
been demonstrated that the introduction of the correlation length for condensate distributions 
produces much smaller values for the first moments of the pion WF compared to the CZ values 
(see also [11], [12]). This leads to a form of the pion WF strongly different from the CZ 
two- hump form and a little wider than the asymptotic form <^"*(x) = Qx{l — x). 

In our previous paper [1] we have obtained directly the forms of the pion and its first 
resonance WFs, using the available smooth ansatze for the correlation functions of the non- 
local condensates and without any suggestion about the quark-hadron duality for the spectral 
density. This program has been suggested and realized in [2], and we have developed alternative 
methods of extracting WFs from this sum rule for a non-diagonal correlator [1]. It should be 
mentioned that we do not need the principle of stability over the Borel parameter M^, and the 
Borel SR appears only at an intermediate stage. 

Here, we test our method using an exactly solvable model as an example (initiated by the 
two-dimensional quantum- mechanical oscillator). We demonstrate the validity of the method 
and also suggest a pure algebraic procedure for extracting the masses and WFs of the particle 
under consideration. We investigate the stabihty of the procedure to perturbations of the 
theoretical part of the sum rule and obtain some results to believe that it works well in real 
QCD SR. In application to the pion case, this results in producing not only the mass and WF 
of the first resonance (tt') but also to provide an estimate of the 7r"-mass. 

^Here iJ(0, z) = Pexp{i (it^A°(t)To) is the Schwinger phase factor required for gauge invariance. 
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Our quantitative results, i.e. the details of the WF shapes and the values of masses, are 
dependent on the certain form of the ansatz used for the non-local condensates. This form may 
be obtained in a future theory of the QCD vacuum. In the absence of such a theory the form 
of the ansatz was chosen as a phenomenological "input" (for details see [1]). From this point 
of view, our calculational scheme may be considered as a suitable framework for connecting the 
hadronic properties, on one hand, with the future theory of the non-perturbative vacuum on 
the other. 



2 Non-diagonal QCD SR for pion wave functions and 
the Method of Integral Transform 

2.1 Non-diagonal QCD SR 

The sum rule for the "axial" WFs (p.^..\x) of the pseudoscalar mesons, based on the non-diagonal 
correlator of the axial and pseudoscalar currents 

^,{x) + (^.,(x)e-^'/^%(^..(x)e-^"/^%...^$(^,x) (1) 

has been suggested in [2] and possesses some remarkable properties. For the WFs on the l.h.s. 
of (1) one has 

{e)M ^ t vm{x){2x - ir dx- {e=\ = i , (^^=°).' = (r=°)." = . . . = o. (2) 

Jo 

On the r.h.s. of (1) there appears a single correlation function /(i/) parametrizing the z"^- 
dependence of the non-local quark condensate: 

(: g-(0)?(z) :) ^ (: g-(0)?(0) :) f e^^/^ f(u) du, (3) 

Here D is the covariant derivative, can be treated as the average virtuality of vacuum quarks. 
Equation (1) demonstrates in the most explicit manner that the distribution (fi-jrix) of quarks 
inside the pion over the longitudinal momentum fraction x on the l.h.s. of the equation is 
directly related to the distribution /(i/) over the virtuality u of the vacuum fields on the r.h.s.. 

It should be noted that the sum rule (1) results from the approximations both in the 
theoretical (for a detailed discussion, see [2]) and the phenomenological parts (see [1]). These 
are 

1. a reduction of three-point correlators to two-point ones; 



^ In deriving these sum rules we can always make a Wick rotation, i.e., we assume that all coordinates are 
Euclidean, z"^ < 0. 
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2. a reduction of the quark-quaik-gluon distribution function to a quark-quark one; 

3. neglecting the contributions of higher non-local condensates, like (qGGq), {qGGGq), . . .; 

4. a representation of the phenomenological spectral density p{s) on l.h.s. of the SR as a 
sum of very narrow resonances 

p{s) = S{s) ■ U^Pt,{x) + ^5{s- Hi) ■ fiipi{x). 

i>l 

In the present work we shall use a specific form of the correlation function /(z^), given by 
f{iy) = N,-exp(-Ayu-a-u); (5) 
Na — — A : — : — iiTi (2;) is the modified Bessel function. 

As was established in our preceding paper [1], this ansatz (A2-ansatz there) conforms with pion 
physics. A similar form for the ansatz also appears in the framework of the instanton model 
for the non-perturbative vacuum ^. Here, the parameter ^ 0.2 GeV^ was extracted from the 
meson QCD SR for the heavy-quark effective theory [2]; the normalization constant Nq and 
the parameter a are fixed by Eq.(4). For the average virtuality of vacuum quarks we take the 
usual QCD sum rule value ~ 0.4 GeV^ [14]. 

As is clear from the structure of the r.h.s. of (1), the form of the function on 
the variable x is given by two humps (one centered at xa = Sj^/M"^, where is the point of 
maximum for the ansatz correlation function /(z^), and the other at xa = 1 — sa/M"^) moving 
as changes. When increases, the humps become narrower, higher and more close to the 
boundary points x — Q ox x — 1. However, ^{^,x) is not yet the pion WF: the larger M^, 
the larger the contamination from higher pscudoscalar states vr', tt", . . .. Due to the properties 
summarized in (2), the corresponding WFs should oscillate (see, e.g.. Fig. 11); therefore the 
function $(-^,a;) does not resemble to (/^^^(x) at not sufficiently low M^. 

At low M^, the pion WF dominates in the total sum (^{^^x) (however, one cannot take 
too low because the operator product expansion fails for < A^). When = 0.4 GeV'^, 
it was observed that ^{^,x) is very close to the asymptotic wave function of the pion (see 
[2]). But there is no strict criterion for selecting the value of to determine the WF of the 
ground state (pion). Our method of extracting the WFs and masses of resonances from (1) 
allows one to avoid this problem. 

2.2 The Method of Integral Transform 

In what follows it is convenient to use a new variable r = 1/ instead of the Borcl parameter 
M^. The r.h.s. of SR (1), i.e., <I>(r, x) is then defined for r G [0, yj]- Let us consider this function 

in the whole complex plane of r and define the integral transformation (IT) P{N,to,uj): 
^we are grateful to A. Dorokhov for informing us about this result (private communication) 
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where the integration is performed along the vertical line C — {c — ioo, c + ioo) with c lying 
to the right with respect to any pole of the integrand, a; > and Tq G (0, p-] . 
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Figure 1: Here C is the contour of integration in the complex plane of r; crosses (x) denote 
the positions (in case A - possible) of the integrand poles. 

To obtain the result of this operator action on the initial SR-representation (here //j = mj) 

"fi^) +^<^i(a;)exp(-T/Xi) = $(T,a;), (7) 

i=l 

consider its action on a simple exponential e~'^^. Evidently, if a; > we can close the contour 
C to the left and get (due to the well-known residue theorem) the contribution from the pole 
at r = To; if /i > we can close the contour to the right and, due to the absence of poles in 
that part of complex r-plane, obtain 0; therefore 



P{N, To, a;) exp(-/xr) = e{uj - /i) ■ exp(-/xro) (^l - 
Then, for a; > 0, we have 



N 



(8) 



i>i 



(jj 



N 



We see that by varying oo one can switch on (switch off) successively more and more resonances 
on the l.h.s. of the SR. So, we have obtained a real step division (smooth at > 0) of the set 
of resonances by the parameter lu in (9), instead of the exponential suppression of the whole set 
of resonances with respect to the inverse Borel parameter r in (7). Accordingly, the hadronic 
characteristics {//j, </?j(x)} are determined by the positions of singularities ( and zeros ) of the 
correlation function /(l/r) in the complex r-plane. 
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Let us assume now that the mass values are known /iq = 0, /xi = Hexp ~ 1.7GeV^, .... 
Then, we can obtain an expression for (fi{x) by employing equation (9) with the parameter 
cu — cui < III that corresponds to the saturation of the ground state (p{x) 

ip{x) = P{N,To,uJi)^T,x) = ^n{uJi,x). (10) 

The final result should not depend on the parameter strongly. However, in a real life we 
don't know the position Hi exactly; therefore we may get contamination from the next state. 
This contamination becomes smaller when N increases due to the power suppression near the 
threshold ni (see l.h.s. of Eq. (9)). One should also take into account that the WF form is 
saturated smoothly in the region cu < /ii and its norm reaches unity at the end of the region. 
The curves corresponding to the function ip{x) = ^^{uJi ~ fJ^i.x) are shown in Fig. 2 for 
different values of N. Their form varies very slightly with = 1, 2, 3 (the same happens 
when the parameter tq decreases from 2.5 to 1.5 GeV^, see the discussion in [1]) and is always 
somewhat narrower than the asymptotic one. Note that formula (10) may be easily generalized 
to a state with a finite width (see [1]). The contribution to the pion WF due to finite width 
effects is not significant numerically. 

We now describe the pure algebraic procedure that allows one to extract from the SR (9) 
both the positions (/Xj) and the WFs {ipi{x)). The procedure will be based on the simple step 
representation appearing on the l.h.s. of (9). Indeed, let us consider the difference 

Di^n{^^, x) = *jv(a;, x) - *jv+i(c^, x). (11) 

We see that the same contributions due to (p{x) are cancelled out and this difference becomes 
non-zero just after cu > /ii. 

Di^Niu, x)=Y.9{uj- fii) ipiix) exp (-tq/x,) ^ (l - I^Y . (12) 

i>l U \ LOJ 

Therefore, this difference is a trigger for the first resonance - the function Di^^iuo^x) must 
have a root at c<j = /ii that fixes the resonance (at fixed x). Analogous trigger-like differences 
could be constructed for higher resonances. By induction, let the difference is trigger-like 
for the n-th resonance (i.e., it has a root at a; = //„) and has the following form 

Dn^]si{io,x)^y26{uj - ^Xi)ipi{x)ex^{-TQHi) —— — — . (13) 

Then the difference D^+i is defined recursively by the recurrence relation 

D^+i^n{u;, x) = ( 1 - — ) Dr,^N{u;, x) - D^^n+i{u;, x). (14) 



For instance, for the second resonance the corresponding difference is of the form 

D2^n(u;,x)^ (l - Di^n{ou,x) - Di^N+i{ou,x), (15) 
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and it yields the value of /X2 provided the value of /xi is known from the previous step. Hence, 
with the help of these trigger-hke differences one can easily determine successively the masses 
Hi. After that it is straightforward to determine the WFs (pi{x) itself, i.e.: 

- $i(/ii,a;))exp(ro/ii) 

M^) = {^^^ 

where the value of the parameter uj e {fii, corresponds to saturation of the first resonance. 
One can try the other formula for (pi{x) following from (12), viz. 

V^iW = — ^iv ■ (17) 



ai (i — H±\ 



The function ipi{x) in (17) should not depend on the value of N in an appreciable manner when 
the saturation of the resonance is reached. 



3 An exactly solvable toy model 

The above formulated procedure is based on the understanding of the initial SR (1) as an 
exact equality. But in the real case the connection between the non-perturbative vacuum 
distribution /(z/) on the r.h.s of SR and the pscudoscalar meson properties reflected on by 
the spectral expansion I]j </^i(x) exp(— /Xjr) on the r.h.s. of SR is only an approximation. 
To understand the reliability of the suggested technique (IT) as an approximated equality, it 
should be useful to check it, first, with a model $(r, a;) having an exact spectral expansion. 
After that we deform $(r, x) in some special way to destroy the exact equality in the spectral 
expansion with the initial WFs and masses, and apply our procedure again. Following this 
way, one may investigate the stability of the procedure to perturbations and its suitability. We 
hope, despite the fact that the exact equality is broken, that the procedure will "work" and 
that the new extracted parameters {/x^, <^^(a;)} will be close to the initial values, if the is not 
too strong. Only in this case our procedure makes sense. 



3.1 Formulating the model 

We should hke to test our approach for some simple exactly solvable model. We select one which 
originates from the 2D harmonic oscillator problem in quantum mechanics (for more details see 
ref.[15]). The Borel transformed Green function M{r,x) of the 2D harmonic oscillator (with 
Wo = 1) in the coordinate representation 

M(r,r)^ ^1^^,(0 1 V-^'= (18) 

k=0 

has for r = the form 

777 

M(t, 0) = ^ , , . (19) 
^ ' ^ 27rsinh(r) ^ ^ 



6 



The spectrum of this model (as one can easily see from the geometric progression summation 
formula) is equidistant with a two-fold step 



Ek = 2k + 1. 



(20) 



We can also write down an explicit analytical form for M(r, r) with r 7^ 0, but there is no need 
for this, since we shall work, for simplicity, with a modified quantity To this end, we define 
the function 

which obviously has the spectral expansion we are interested in, viz. 



k=0 



k-2Tk 



(22) 



This toy model has a nice feature: all its resonance WFs are defined through one WF, namely 
the WF of the first resonance ipi{x) — (f){x) which can be selected at free choice. For our 
convenience, we take the following form 



(t){x) = 4a;(l - x), 

i.e., the asymptotic WF of the pion (albeit with wrong normalization). 



(23) 



3.2 IT-technique for the toy model 

For evaluating $Ar(cj, x) in our toy model, we use the well-known residues-technique. Separating 
residues (sec Fig. 1) from the To-pole (A^-q^jv) and the poles Tn{x) = \{x) + iim with \{x) = 
|log(0(a;)) (A$<l>Ar), we obtain 



^n{uj,x) ■ 



exp(— u;ro) 



UJ 



N 



exp(a;ro) 



1 - Lp{x) exp(-2ro) 



A,$. = -g(a;)r(7V+l) ^"P^-"'"°:^(")]> x 



X 



+ E 



cos [i^nuo + (A^ + 1)0„] 

{N+l)/2 



2(ro-A(x))^+^ nVi[(ro-A(a;))^ + 7r 



2^2 



Here 0„ = arctan [vrn/ (tq — A(a;))]. 

The l.h.s. of SR (9) in this model has the form 



1 + ^ e (u; - 2A;) (t>{xf 1 exp (-2A;to) 

fe>i V ^ / 



By this way we obtain the exact SR in the form 

(24) = (27). 



(24) 
(25) 



(26) 



(27) 
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Saying "exact" we mean, that this SR is vahd for all x e [0, 1] and for all uj. Moreover, it 
holds also for all choices of the initial WF (f){x) ! This behavior is shown exphcitly in Fig. 3, 
where the step-like structure on the l.h.s. of SR (see (27)) is reproduced by its r.h.s. (see 
(24)). Our trigger- like differences exactly determine masses /i^ — 2A; and the corresponding 
WFs (pk{x) = 4>{x)''. 

3.3 Deformations of $(r, x) and the IT-procedure 

All formulas of the preceding subsection are exact (they are in some sense identities). But in 
real QCD problems we have no such regular behavior of the theoretical part of the SR and one 
can use the spectral expansion only approximately. Moreover, we cannot be sure about the 
adequacy of ansatz (5) for the true non-perturbative amplitude in the whole region of variables. 
Having this in mind, let us inspect the stability of the suggested procedure under perturbations 
of the r.h.s. of (9). 

The first problem here is the following: how can one deform the original function $(t, x) 
while retaining the possibihty of controlling the degree of this deformation? 

We apply the following method. As we know, the function $(t, x) can be expanded in a 
series of poles over the variable r: 

1 1 °° 1 

*<^' « " l-^Wexp(-2.) = 2 + W^^y '''' 

If we exclude in this expansion a set of poles, we obtain some deformed function. It should be 
stressed that this function does not have the same spectral expansion as the original one. Let 
us consider our IT-procedure with these deformed functions. The exact result for the difference 
Di^i is reproduced in Fig. 4. The deformations we have determined will destroy the exact 
step-structure on the l.h.s. of (9) and (13). This procedure proceeds as follows. 

1. At first we exclude all the poles t„(x) with \n\ > K (the corresponding deformed $ will 
be denoted as ^{u!,x; K)). Fig. 5 contains the curve for Di^i{uj , x; K = 6). As we can 
conclude in this case, the situation is fine: the main deviations from 6{uj — 2) are located 
near the origin (a; ~ 0.2). The locations of resonances are well-determined. 

2. At the next step we exclude also the last six poles (|n| = 3, 4 and 5). The corresponding 
curve for Di^i{u,x; K = 3) is reproduced in Fig. 6. The situation is still good: though 
deviations from 9{u! — 2) are larger relative to K = 6, they do not extend over cu ~ 0.2. 
The locations of resonances are also well-determined and our "triggering" is still vahd. 

3. The last step is to exclude the same poles as in the first case, but in addition also those 
with n = ±1. The curve for Di^i{u!,x; K = 6 -|- 1) is reproduced in Fig. 7. Now the 
pattern is changed: deviations form the initial form are large and they reach u — 1. 
Indeed, this new form is rather close to ^ (a; — 1): there exists the zero for a; = 1 where the 
original function has no resonances at all. Going through the procedure with this value 
(Hi = 1) for the mass of the first resonance, we obtain the new spectrum fik = k with 
the new resonance WFs {(^^(x)} 7^ ]0*^(a;) [. The crucial indication for the correctness of 
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such an interpretation is the smallness of the deficit of the initial SR 

AsRmT,x;K ^6 + 1)} = ^T,x;K = Q + l) (29) 

- [(po{x) + (pi{x) exp(-/xiT) + ip2{x) exp(-/X2T)] + . . . . 

We evaluate this quantity for the new and for the initial spectra and obtain a clear signal 
that the new one provides the genuine description of the SR (see Fig. 8 and 9). 

As to the extracted wave functions ipk (x) , we can characterize their evolution under deforma- 
tions by the following statements: 

• opposed to the exact solution, {ipk{x)} become a;-dependent when cu e (/Xj; yUj+i]. They 
tend to the exact ones with cu more close to the right end of the corresponding interval 
(see Fig. 10); 

• the difference between (pk{x) and the exact WF of the A;-th resonance becomes stronger 
for higher k (i.e. for highly exited states; see Fig. 10). 

To conclude these investigations: the method really works, even for strong deformations of 
the structure of the initial spectrum as the last example of deformation demonstrates. Moreover, 
we see that deformations of the analytical structure of the initial $(r, x;0) in the complex r- 
planc far from (near) the main pole r = \{x) produce small (large) deformations of the spectrum 
and the WFs. Therefore, the knowledge of the analytic properties of the function $(r, x) in a 
region in the complex r-plane when includes the pole r = X{x) and a few others appears to be 
sufficient for obtaining information on the first few resonances and their WFs. We hope that 
an analogous situation appears in the SR for pion. 

4 IT-procedure for the pion and its resonances 

Applying the algebraic procedures to the IT of SR (1) we arrive at the following results: 

1. For uj — — 1.8 GeV^, the difference Di^i{ll!.^,, 1/2) is equal to zero and with increasing u 
it becomes negative (this means that (/9i(l/2) < 0). This value of the first-resonance mass 
is quite reasonable for the pion case, the experimental one being /i^^^ = m^, ~ 1.7 GeV 
(at the full width F^r' = 0.2 0.6 GeV). The pion WF (p{x), corresponding to this uj^ has 
been obtained from expression (10) in subsection 2.2 (see Fig. fig-A2-123). Its width is 
traditionally characterized by the value of the second moment - (^^)7r ~ 0.17. Another 
characteristic often appearing in form factor calculations is the integral Jq (p{x) / xdx ^ 2.8. 
Note that the method does not work properly in the neighborhood of the cndpoints x = 
and a; = 1 of the WF. Therefore, to estimate the above values, one should exclude these 
regions in the integral. 

It should be emphasized that ^i{u,x) imitates the l.h.s. of the SR as a function of u; 
rather well, despite its different origin. Note here that the position of the root a;* depends 

also on the value of x: in the region 0.35 < a; < 0.65 this dependence is rather weak. We 
don't see the step structure on the l.h.s. of Eq. (12) - the structure is destroyed due to 
the employment of too crude approximations in the theoretical part of the SR. Following 
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our experience in the "exactly solvable toy model" considered in subsection 3.3, we may 
hope, nevertheless, that the roots of the few trigger differences -Di,2,...^Ar(<^*, 1/2) yield 
reliable approximations for the masses. 

2. The saturation by the first resonance is reached near u; 3 GeV^. The curves correspon- 
ding to the resonance (pi{x) (for u; 3 GeV^) are shown in Fig. 11. 

3. The last information we can extract in this case is the position of the second resonance: 
the equation D2^i{uj**,x) — gives us the range of values a;** — IJ,2 — 3-7- 4GeV^, which 
is converted to the mass value m2 — 1.86 ± 0.14GeV. The experimental parameters for 
this resonance are still in question in the Particle Data booklet (April 1994). However, 
the recent measurements in [17] certainly provide for the mass m^r" ~ 1.78 ± 0.007 GeV 
(at full width Ftt" = 0.16 GeV), which is rather close to our estimate. 

5 Conclusion 

We have considered a model SR for the pion and the pseudoscalar resonance WFs based on the 
non-diagonal correlator introduced in [2]. The theoretical side (r.h.s.) of this sum rule depends 
only on the non-local condensate. We have tested the approach proposed in our recent paper 
[1], which enables one to extract WFs and masses from this SR, using the exactly solvable 
toy model as an example. We demonstrated the vahdity of the method and suggested a pure 
algebraic procedure for extracting the masses and WFs relating to the case under investigation. 
We investigated also the stability of the procedure under perturbations of the theoretical part 
of the sum rule. We obtained as a result that the most crucial domain in the complex plane of 
T = is the neighborhood of the poles lying on the real axes - the perturbations near these 
points essentially reorganize the spectral expansion under investigation. Applying this method 
to the pion case and using one of the ansatze given in [1] not only the mass and WF of the 
first resonance (tt') have been derived, but also the mass of n". Our results confirm the main 
conclusions about the shapes of the WF of the pion and its resonances, obtained in [2]. 

It should be mentioned that the specific form of the ansatz for the correlation function 
/(z/) plays an important role in determining concrete values of masses and shapes of the WFs 
corresponding exactly to the pion resonances (on the contrary, the WF of the pion ground 
state is not too sensitive to the shape of the ansatz). We have chosen our ansatz (5) among 
others (see e.g. [1]) just due to this reason. So, one may consider it as a first step towards 
solving the inverse problem, namely to obtain the vacuum condensate properties from hadron 
phenomenology. 
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Figure 2: Forms of the WF of pion (pnix) extracted from ~ 1-5, x) for — l(line), 

N = 2 (dashed) and N = 3 (dotted) 



Figure 3: Step-hke structure of $i(a;, x) with x — fixed = 0.5 for the toy model case. 
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Figure 4: Difference Di^i{u},x) with x — fixed = 0.5 for the exact r.h.s. of SR for the toy 
model case. 



Figure 5: Difference Di^i{uj,x, K — 6) with x — fixed = 0.5 for the deformed {K — 6) r.h.s. 
of the SR for the toy model case. 
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Figure 6: Difference Di^i{uj,x, K — 3) with x — fixed = 0.5 for tfie deformed {K — 3) r.li.s. 
of tlie SR in tlie toy modei case. 



Figure 7: Difference Di$i(cj, x, K = 6 + 1) witfi x — fixed = 0.5 for tfie deformed {K = 6 + 1) 
r.li.s. of tlie SR in the toy model case. 
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Figure 8: The deficit of tlie initiaf SR Asr{^{t, x; K — 6 + 1)} for tlie new spectrum /ik — k 
in tlie toy model case. 



Figure 9: The deficit of the initial SR Asr{^{t, x;K — 6 + 1)} for the old spectrum fj,k — 2/c 
for the toy model case. 
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Figure 10: The difference between ipi{x) and the exact WF of the i-th resonance for the i — 1,2 
for the deformed {K = 6) toy model case. 



Figure 11: Form of the WF of the first resonance (p-n-'{x) extracted from ^ 3 GeV^, x) 

for the pion case. 
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